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The measurement of the beauty quark production cross section in deep 
inelastic scattering with the ZEUS detector at HERA is presented. The 
cross section for beauty quarks, measured in the kinematic range > 2 
GeV^ and 0.05 < y < 0.7, and for events with at least one muon and 
one jet in the Breit frame, is somewhat above the next-to- leading order 
(NLO) QCD prediction, but in agreement within the experimental and 
theoretical uncertainties. 

1 Introduction 

Deep inelastic scattering (DIS) offers the unique opportunity to study the pro- 
duction mechanism of heavy quarks via the strong interaction in a particularly 
clean environment where a point-like projectile, a virtual photon with a vir- 
tuality Q^, collides with a proton. A first measurement of the total inclusive 
6-quark cross section in the DIS region has been released by the HI collabora- 
tion [1]. Recently, the ZEUS collaboration has measured the production cross 
section of 6-quarks in DIS at HERA, in the reaction with at least one hard jet 
in the Breit frame^ and a muon in the final state (ep —> ebbX e jet ^ X) [2]. 

Due to the large 6-quark mass, muons from semi-leptonic b decays usually 
have high values of p^'' , which is the transverse momentum of the muon with 
respect to the axis of the closest jet. This allows a statistical separation of the 
signal from the background. 

*Now at the Physik Institut der Universitat Ziirich-Irchel, 8057 Zurich, Switzerland. 

^In the Breit frame, defined by 7 -|- 2xP = 0, where 7 is the momentum of the exchanged 
photon, X is the Bjorken scahng variable and P is the proton momentum, a purely space-like 
photon and a proton collide head-on. 
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2 Event selection and extraction of the b decays 

The data used in the measurement were collected during the 1999-2000 HERA 
running period, corresponding to a total integrated luminosity of 60 pb~^ and 
at a center-of-mass energy of 318 GeV. A detailed description of the ZEUS 
detector can be found elsewhere [3]. 

The event sample was selected in four steps: 

• DIS events were selected by requiring a well reconstructed outgoing po- 
sitron with energy greater than 10 GeV, > 2 GeV^ and inelasticity 
0.05 <y < 0.7, where y = Q'^/xs. 

• Muons were identified by requiring a reconstructed segment both in the 
inner and outer part of barrel and rear muon chambers, where the rear 
direction is defined by the incoming positron. The reconstructed muons 
are matched in space and momentum with a track found in the central 
tracker. In addition, a cut on the muon momentum > 2 GeV is applied 
together with a cut on the muon polar angle 6^, 30° < 6^^ < 160°. 

• Hadronic final-state objects, reconstructed from tracks and energy de- 
posit in the calorimeter, were boosted to the Broit frame and clustered 
into jets using the fcy clustering algorithm [4]. Selected events had at 
least one jet with E^"^^^* > 6 GeV and within the detector acceptance 
— 2 < rj^"-^ < 2.5, where r?^"'' is the pseudorapidity in the laboratory 
frame^ . 

• The muons in the sample were associated with a jet using the kr algo- 
rithm information, where the associated jet was not necessarily the jet 

satisfying the cuts above. To ensure a good reconstruction of the asso- 
ciated jet, it was required to have transverse energy in the Breit frame 

^Breit ^ 4 QgV. 

After all selection cuts the final data sample contains 836 events. 

To correct the results for detector effects and to extract the fraction of 
e\(>utri from /) dccaxs llio RAPGAP Monte Carlo program [5] was used. The 

"^Thc pseudorapidity r; is defined as r; = — ln(tan$/2), where is the polar angle measured 
in with respect to the proton beam direction 



3 



light-flavour and charm-quark RAPGAP samples were mixed according to the 
relative luminosities and the 6-quark sample was added according to the beauty 
fraction determined from the analysis of the p^*^' distribution. To determine 
the fraction of events from b decays in the data, the contribution from light- 
plus-charm flavours and beauty were allowed to vary, and the best mixture was 
extracted using a binned maximum-likelihood method. The measured fraction 
of b decays is (25 ±5)%, where the error is statistical. This mixture describes 
well the shape of all event kinematic variables. 

3 Results 

The total visible cross section, fTj,^, was determined in the kinematic range 
Q2 > 2 GcV^, 0.05 < y < 0.7 with a muon with > 2 GeV and 30° < 
6*^ < 160° and at least one jet in the Breit frame with i;^""*'* > 6 GeV and 
-2 < 77-^°^ < 2.5. The cross section is 

a^lie+p e+bbX e+ jet X) = 38.7 ± 7.7 (stat.)+5;o(syst.) pb. (1) 

The NLO QCD predictions were evaluated using the HVQDIS program [6]. Frag- 
mentation of 6-quarks into hadrons was performed using the Peterson function 
[7] with e = 0.002 as suggested in Ref. [8]. The semi-leptonic decay of 6-hadrons 
was modeled using a parameterisation of the muon momentum spectrum ex- 
tracted from RAPGAP. The 6-quark mass was set to nt), = 4.75 GeV and the 
renormalisation and factorisation scales to /i = ^JQ"^ + Ami. ^^"^ CTEQ5F4 
proton parton densities [9] were used. The NLO QCD prediction is 2813:5 pb, 
where the error was estimated by varying the scale /i by a factor of 2 and the 
mass nib between 4.5 and 5.0 GeV and adding the respective contributions in 
quadrature. A more detailed discussion of the QCD prediction and the related 
uncertainties for this process can be found in Ref. [10]. The measured total 
cross section is somewhat above but agrees with the NLO prediction within the 
experimental and theoretical uncertainties. The program CASCADE [11], based 
on the CCFM evolution equations [12] and on gluon densities depending on the 
transverse parton momenta emitted along the cascade, gives a prediction of 35 
pb, which is in good agreement with the data. 

The difl[erential cross sections were calculated in the same kinematic range 
as the total cross section by repeating the fit of the pl^^ distribution in each 
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bin. Figure 1 shows the differential cross section as functions of and x 
compared to the NLO. The NLO predictions agree reasonably well with the 
data. Figure 2 shows the same differential cross sections compared with the LO 
Monte Carlo programs. The LO QCD prediction folded with the DGLAP QCD 
[13] evolution equations (RAPGAP) underestimates the measured cross section. 
CASCADE is in good agreement with the data. 

ZEUS ZEUS 
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Figure 1: Differential beauty cross section as a function of (left) and as a 
function of x (right), compared to the NLO QCD calculations. The shaded 
bands show the uncertainty of the theoretical prediction {see text). 



4 Conclusions 

The production of beauty quarks was measured in the process e~^p — ^ e"^ + + 
jet + X with the ZEUS detector at HERA. The NLO QCD predictions agree 
with the measured total cross section and differential cross sections as func- 
tions of and the Bjorken variable x within the experimental and theoretical 
uncertainties. The CASCADE Monte Carlo program gives a good description of 
the measured cross sections. 
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Figure 2: Differential beauty cross section as a function of (left) and as a 
function of x (right), compared to the LO QCD Monte Carlo programs. 
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